Singularities and Avalanches in Interface Growth with Quenched Disorder 
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A simple model for an interface moving in a disordered medium is presented. The model exhibits 
a transition between the two universality classes of interface growth phenomena. Using this model, it 
is shown that the application of constraints to the local slopes of the interface produces avalanches 
of growth, that become relevant in the vicinity of the depinning transition. The study of these 
avalanches reveals a singular behavior that explains a recently observed singularity in the equation 
of motion of the interface. 
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The problem of interface motion in disordered media 
has attracted considerable attention recently. By means 
of numerical, analytical, and experimental studies EJ it 
has been observed that scaling theory can be used as 
the framework to understand the dynamical properties 
of the interface. In a typical realization of the problem 
a (i-dimensional interface moves in a (d+ l)-dimensional 
disordered medium, driven by an external force F . The 
interface is assumed to be oriented along the longitu- 
dinal a;-direction, and is specified by the height y(x,t). 
For small forces the interface is pinned by the random 
quenched impurities of the medium, and the average ve- 
locity is zero. Above a critical value F c , the external force 
overcomes the effect of the impurities, and the interface 
moves with a finite velocity. Near the threshold, the ve- 
locity scales as v ~ f e , where / = F/F c — 1 is the reduced 
force, and 9 the velocity exponent. 

The roughness exponent a is defined by the relation 
W(L) ~ L a , where W is the width of the interface and 
L is the system size. The exponent a characterizes the 
different universality classes of interface growth phenom- 
ena. The universality classes can also be classified [|] ac- 
cording to the behavior of the coefficient A of a nonlinear 
term of the form A(Vy) 2 in the equation of motion of the 
interface. One of the universality classes is described by 
an equation of the Kardar-Parisi-Zhang (KPZ) type || 
with a quenched random field rj(x,y) 



d t y(x, t) = W 2 y + \{Vy) 2 + rj(x, y) + F, 



(1) 



with the coefficient A diverging at the depinning transi- 
tion as 



a(/) ~r , 



(2) 



where (f> ~ 0.64 0. The roughness exponent at the depin- 
ning transition in (1 + I)-dimensions is a — 0.63, and can 
be obtained by a mapping to directed percolation 
We refer to this universality class as directed percolation 
depinning (DPD). 

A second universality class is described, at the depin- 
ning transition, by an equation of the Edwards- Wilkinson 
tyP e @ with quenched disorder 



Models with A = (for any force), or A — > (when 
F — > , corresponding to a negative exponent <f> in Eq. 
(H)), belong to the universality class of Eq. (||). Analyti- 
cal and numerical [^|j9| studies in (I + I)-dimensions 
yield a roughness exponent a ~ 1.23 — 1.25. This univer- 
sality class is referred to as quenched Edwards- Wilkinson 
(QEW). 

In this paper we present a simple model for an elas- 
tic interface that exhibits a transition between the two 
universality classes of interface growth in a disordered 
medium. The DPD result is obtained when lateral or 
longitudinal fluctuations of the interface height are al- 
lowed and the growth rule is restricted by a generalized 
solid-on-solid (SOS) condition. On the other hand, when 
the constraint is absent, the QEW universality class is 
obtained. The generalized SOS condition is a universal 
feature appearing in all the models of the DPD univer- 
sality class. It is analogous to the restricted SOS condi- 
tion (\y(x ± 1, t) — y(x, t)\ < 1) presents in simple models 
of growth with time- dependent noise 0. In these type 
of models, the constraint is commonly associated with 
lateral propagation of growth that generates nonlinear- 
ities in the equation of motion. Moreover, it is shown 
that, in the presence of quenched disorder, the general- 
ized SOS condition also generates avalanches of growth. 
These avalanches are irreversible growth events that be- 
come relevant near the depinning transition. We show 
that they present a singular behavior that is responsible 
for the divergence in the coefficient A. 

Consider a one-dimensional elastic interface moving 
in a two-dimensional disordered medium of lateral size 
L. A discrete model for such an interface is defined 
in the square lattice by the height values {y k }k=i,--,L- 
The interface moves under the influence of an external 
force F — (0,F) and fluctuates in the longitudinal or re- 
direction, and transversal or y-direction. The strength of 
these fluctuations is controlled by the elastic constants v x 
and v y , respectively. For the A;-th column and at a given 
time, the velocity v = (v%, v y ) is calculated according to 



v v (Vk+i + Vk-i - 2^fe) + r](k, y k ) + F 
-2v x + rj(k,y k ). 



(4) 



d t y{x,t)=y 2 y + n(x,y)+F. 



(3) 



The term rj{k, y k ) represents the strength of the random 
force, and mimics the effect of quenched disorder. It is 



1 



a random field with a uniform probability distribution 
between [—5, 8], and represents a repulsive force for posi- 
tive values and an attractive or pinning force for negative 
values. A site at the interface moves forward or laterally 
only when the respective velocity is positive. Therefore, 
the k-th column and its nearest neighbors are updated 
in the following way 



Vk 

Vk+i 
Vk-i 



Vk + 1 if 
Vk if v, 

Vk if v. 



Vy > 



> 

> 0. 



(5) 



The first equations in (Q) and (||), correspond to a 
model introduced to study an elastic interface described 
by the QEW equation || , and they are equivalent to the 
discretization of Eq. (|^) . The last equations in (0) and (||) 
are the simplest generalization of the model to include 
longitudinal motions [fll|| . A lateral motion to a nearest 
neighbor column is allowed if the neighboring column is 
smaller than the column considered. After the lateral mo- 
tion, the interface becomes a multi- valued function of the 
longitudinal coordinate x (see Fig.|l|(a)). The generalized 
SOS condition is then applied in order to transform the 
interface into a single- valued function, by defining the in- 
terface with the highest value of the height. This growth 
process can be thought as a transversal avalanche of 
growth in the neighboring column, as shown in Fig^(6). 
This growth event occurs regardless the value of the noise 
in the neighboring column. 

We perform numerical simulations on a square lat- 
tice with L = 1024. Helical boundary conditions tjl = 
yi + m L, where m = (Vy) is the average external tilt, 
are imposed on the interface in order to study the in- 
terface velocity as a function of the tilt Figure || 
shows the tilt-dependence of the velocity of the interface 
for the value v x = 0.2. In the following, we use the values 
v y = 1.0 and 5 = 3.0 since our results do not depend 
of these parameters. We see that the parabolas become 
steeper when F — > F c + as in Eq. (||), indicating that the 
model belongs to the DPD universality class. 

The model presents a transition at a critical value v Xc . 
If the value of v x is increased such that v x > v Xc = 5/2, 
then the QEW result is found. Specifically, we find that 
v is independent of m, indicating that A = for any 
force. The fact that for v x > 5/2 the ^-component of the 
velocity v x is always negative explains this transition: 
longitudinal motions cannot occur and one recovers the 
model of Ref. [0 with the corresponding behavior of A. 

The different parabolas obtained for a given value of 
v x < v Xa , can be rescaled using the scaling ansatz M 



v(m,f)~f° 9{m 2 /f+% 



(6) 



where g is a universal scaling function. Support for (g) 
is provided by the data collapse shown in Fig. ^, where 
we replot the data of Fig. || and the data obtained for 
v x = 0.6, v x = 1.0, and v x = 1.4 |l4| ]. The scaling func- 
tion g becomes flatter as v x — > v~ , indicating that the 



prefactor of A goes to zero, even though, the singularity 
associated to / is still present as long as v x < v Xa . 

In the following, we argue that the divergence in A 
is explained by the singular behavior of the size of the 
avalanches of growth. 

According to (^|), the interface can advance in two in- 
dependent ways. One way is via a transversal motion (if 
v y > 0), and the other is via a longitudinal motion (if 
v x > 0) plus the transversal avalanche in the neighboring 
column. In order to determine the relevant growth mech- 
anism near the depinning transition, we study the mean 
value of the number of longitudinal and transversal mo- 
tions per unit time n x and n y , respectively, as a function 
of the force. For v x — 0.2, we find 



n-x(f) 
%(/) 



(7) 



with r) x ~ 0.60 and -f y ~ 0.78. Both quantities go to zero 
at F c , since the velocity vanishes at the depinning tran- 
sition. However, the ratio n x /n y ~ /~ 018 diverges at F c , 
indicating the relevance of longitudinal fluctuations for 
the motion of the interface at the depinning transition. 
This fact can be explained as follows. For forces close to 
F c the velocity is almost zero and the interface moves in 
a very irregular way, jumping from one metastable pin- 
ning configuration to another. In this "jerky" motion the 
interface takes advantage of longitudinal motions, rather 
than transversal ones, to surround and overcome the im- 
purities. Thus, avalanche events become relevant for the 
motion of the interface only in the vicinity of the depin- 
ning transition. 

We also study the mean value of the size of the 
avalanches produced by the generalized SOS condition 
per unit time, (s), as a function of the tilt and for dif- 
ferent forces. Figure |] shows the results. As it turns out, 
(s), as well as the velocity, has a parabolic dependence 
on m. These parabolas become steeper as F — » F+, and 
we can fit (s) to 



(s) = s + X s m 2 



(8) 



with X s ~ /~^, and tfi ~ 0.64 the same exponent as in Eq. 
(||) . The parabolas can also be rescaled using the scaling 
ansatz of Eq. (^) with the same value of 9 used for the 
velocity curves. 

As shown in Fig. [l](c) , the size of an avalanche is larger 
for the tilted interface than for the untilted one. This ex- 
plains the increase of (s) with the tilt, exemplified in Fig. 
([|) for a given fixed force. Moreover, since the relative oc- 
currence of lateral motions and avalanches is larger near 
the depinning transition than away from it, the same ex- 
ternal tilt will cause a larger increase in the average (s) 
near F c than for F ^S> F c . Thus, the coefficient A s , that 
measures the variation experienced by (s) due to a change 
in the average tilt of the interface, increases its value as 
F — > F c and the parabolas become steeper. 

Notice that the relevance of longitudinal motions and 
avalanches of growth near F c implies that the velocity 
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is determined by the size of the avalanches, v oc (s), so 
that A oc A s . Thus, the singularity of the coefficient A can 
be explained by the same divergence observed in A s . The 
motion of the interface near the depinning transition is 
entirely dominated by the avalanches of growth produced 
by the generalized SOS condition. 

The same behavior can be predicted for the other mod- 
els of the DPD universality class. All these models share 
a constraint in the growth rule of the interface height 
that generates avalanches analogous to the ones in our 
model. In the model of Ref. a slope constraint is ap- 
plied to the interface that implies a readjustment of the 
height regardless of the value of the noise. The so-called 
erosion of overhangs in the model of Ref. H corresponds 
to our generalized SOS condition. And the model of Ref. 
flBf presents a restricted SOS condition that produces 
avalanches of growth (see Ref. [|l6| ) . 

Other variants of the model (Q) are also studied. First, 
we study a model that includes lateral motions but re- 
moves the generalized SOS condition. The model can be 
better understood as a set of L-beads fluctuating in both 
directions and interacting through elastic springs |p7| . 
We calculate the velocity for this model as a function of 
the average tilt and for different forces. We find a small 
value of A for F 3> F c and A — > when F — ► F+. There- 
fore, the model falls in the QEW universality class. The 
nonlinearity observed in this model is of kinematic origin 
(A oc v) ||, in contrast to the behavior of A in the DPD 
universality class, for which A is coupled to the external 
force and diverges at F c . These results suggest that lat- 
eral motions alone generate nonlinearities only when the 
motion of the interface is so fast that the noise can be 
regarded as a time-dependent noise . At the depinning 
transition, where the quenched disorder is relevant, the 
QEW result is recovered. These results confirm that the 
generalized SOS condition is the origin of the singularity 
in A. 

We also study other mechanisms, which are believed 
to generate nonlinearities in the equation of motion [ p"8[ : 
(i) an external fixed force acting in the normal direction 
of the interface, analogous to the Lorentz force found in 
the motion of flux lines in type II-superconductors Jl9| ; 
and (ii) an anisotropic random medium characterized by 
an anisotropic random force H]. In order to check this 
hypothesis, we simulate the model that allows for longitu- 
dinal motions but without the generalized SOS condition 
pof . In addition, we generalize the definition of the veloc- 
ity vector by including: (i) an external force, with fixed 
magnitude F, acting in the direction of the local normal 
vector h of the interface: F = (F x ,F y ) — Fn; and (ii) 
we model the effect of the anisotropy in the medium by 
introducing a noise vector if = (rj x ,riy), where r\ x and rj y 
are independent random fields with amplitudes S x and 
S y respectively (S x ^ 5 y for an anisotropic medium). We 
find that the model does not belong to the DPD univer- 
sality class. Under the limitations of our model, we find 
that the anisotropy of the medium does not generate the 



divergence in A. 

In conclusion we present a simple model for an elas- 
tic interface moving in a quenched disordered medium, 
which captures the relevant features of the two univer- 
sality classes. The origin of the singular behavior in the 
equation of motion of the DPD universality class is ex- 
plained by a simple microscopic constraint imposed to 
the growth of the interface that generates irreversible 
avalanches of growth. 
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FIG. 1. Sketch of the effect of lateral motions of a site at 
the interface: (a) -(b) for an untilted interface, and (c) for a 
tilted interface with m = 1 as defined in the text. A lateral 
motion (a) produces an effective avalanche (b) in the nearest 
neighbor column due to the generalized SOS condition. The 
size of the avalanche, s, is larger for the tilted interface than 
for the untilted one, as can be seen comparing (b) and (c). 

FIG. 2. Plot of the velocity as a function of the average tilt 
of the interface for values of the reduced force from / = 0.19 
(bottom) to / = 1.38 (top). The parameters are v x = 0.2, 
v y = 1.0 and 8 = 3.0. Results are averaged over 70 indepen- 
dent realizations of the disorder. The "closing" of the parabo- 
las shows that a diverging nonlinear coefficient is present in 
the equation of motion. 

FIG. 3. Data collapse of the data of Fig. 2 (v x = 0.2) and 
the data corresponding to u x = 0.6, v x = 1.0 and v x = 1.4 
(shown from bottom to top, respectively), plotted according 
to the scaling relation of Eq.(6). The scaling function becomes 
flatter as v x — > u Xc = 5/2, indicating the transition to the 
QEW universality class. Each set of curves is shifted for clar- 
ity. 

FIG. 4. Plot of the average size of the avalanches produced 
by the generalized SOS condition as a function of the tilt, for 
the same forces and parameters as in Fig. 2. 
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